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Abstract 

In a first part we study the <^ p+1 -field theory from the classical point of view. Using 
Butcher series we compute explicitly the perturbative expansion of the solutions and we prove 
that this expansion converges if the coupling constant is small enough. Then we show that we 
can formally recover the Heisenberg's interacting quantum field directly from this expansion. 
In other words we write the Heisenberg interacting field as a Butcher series. 
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Introduction 

Butcher series are sums indexed by planar trees introduced by J. C. Butcher [6] in order to study 
and classify [8] Runge Kutta methods in numerical analysis. But they appear naturally in a large 
class of problem (see e.g. [9], [4]) including quantum field theory (see e.g. the paper of Ch. Brouder 
and A. Frabetti [2], [3], [5] about Q.E.D). Moreover Ch. Brouder noticed that the structure which 
underlies the original Butcher's calculations [6] is exactly the Hopf algebra of rooted trees defined 
by D. Kreimer in his paper about renormalization [11] and he shows that finally renormalization 
can be seen as a Runge-Kutta method. Hence it seems that the work of J.C. Butcher and perturba- 
tive quantum field theory are closely related. In this paper we show that formally the Heisenberg's 
interacting quantum field can be written as a Butcher series. Hence we link the perturbative clas- 
sical field theory with the perturbative quantum field theory. 



Butcher series provide a precise description of the solutions of a non linear problem. In this 
paper we focus on p+1 -theory but our results can be generalized to any theory. First, we explain 
how Butcher series give the explicit perturbative expansion of the classical field and we prove that 
this expansion converges if the coupling constant is small enough. Then, we formally quantize the 
expansion and we show that we recover the Heisenberg picture of the interacting quantum field. 

Let p be a non negative integer such that p > 2. Then consider the equation over scalar fields 
ip : R n — ► M 

(□ + m 2 )ip + Xip p = (K-Gp) 

where m > is the mass, A is a real parameter (the coupling constant) and □ denotes the operator 
d ( ® U yi — J2j Q^jyi ■ We show that the solution of (K-Gp) writes as a power series indexed by planar 
trees. 



Planar trees are rooted trees drawn into the plane with the root on the ground. The external 
vertices are called leaves and the other are called internal vertices. We denote by \b\ the number 
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of internal vertices of a planar tree b. There is a unique planar tree with no internal vertices: this 
is the planar tree reduced to a root, let's denote by o this trees. We say that a planar tree is a 
p-tree if and only if each internal vertex has exactly p childrcns. Let us denote by T(p) the set of 
p-trees. Finally given (61, . . . , b p ) a p-uplet of p-trees, we can define another p-tree B + (bi . . . b p ) 
by linking a new root to the roots of 61, . . . , b p . 

Given the Cauchy datas {ip°, tp l ) G H q+1 (R d ) x H q (R d ), we show (theorem 2.1) that the solution 
tp € C 1 ([0,T},Hi(R d ))r\C 2 ([0,T},Hi- 1 (R d )) of (K-Gp) such that ^(0,«) = if and §f(0,») = if 1 
is given by the power series 

b£T(p) 

where the coefficients (0(fe)) beT ( p ) are such that 

{(j)(o) is the solution of (□ + m 2 )0(o) = with Cauchy datas ((p°, ip 1 ) 
4>(B + (bi . . . bp)) is the solution of (□ + m 2 )(p = —(p(bi) ■ ■ ■ 4>(b p ) with zero Cauchy datas 

The power series (B ) converges in the C x ([0, T],H q (R d )) n C 2 ([0, T],H q ~ 1 (R d )) topology if A is 
small enough, it is called the Butcher series. 

Remark 

The function </>(o) is the (classical) free field corresponding to the Cauchy datas (ip°, ip 1 ) at time 
t = 0. Wc can reformulate the definition of <t>(B + (b\ . . . b p j) by 

<A(B+(6 1 . . . b p ))(x) :=- [ dyG ret (x - y)^b x ){y) ■ ■ ■ 4>{b p )(y) 
Jp+ 

where G ret (z) denotes the retarded Green function of the linear Klein-Gordon operator (□ +m 2 ) 
and where P+ is the half space P+ := {(t, ~x) £ R x R d \ t > 0}. 

Using the last remark we can formally quantize the Butcher series. We use the notation of the 
book of E. Pcskin and Daniel V. Schroeder [12] p. 83. i.e. we consider the operator <pi(x) acting 
on Fock space defined by 

; (2ir) d J Rd ^%T k \ k J x°=t-t 

Here k ■ x denotes the quantity k ■ x := io k t — (k, ~~x) where (k, ~~x) denotes the usual scalar product 
of R d and u>k the quantity ujk := {\k\ 2 +m 2 ) 1 / 2 . The operator 4>i{x) is the free field on space-time. 

Then by analogy we define formally the family of operator (<f>(b))beT(p) by 

\ HB+(h . . . b p ))(x) := - d d+1 yG ret (x - y)^(h)(y) ■ ■ ■ $(b p )(y) 

Note that this definition is formal since operator product is ill defined on the diagonal. 

Using the commutation rules of the 0/(x)'s, we get the following theorem: 
Theorem 3. 1 for all t > to and ~x <E R d we have formally 

(p{t,lc) := Y A |b| 0(6)(x) = U\t) <f>i(t,lc) U(t) 

b£T(p) 



where U(t) denotes the operator 



ft rTi rT, 

U(t):=J2(-i*) a dn / dr 2 ... 

q >q Jto J t J to 



dT a H I (T 1 )H I (T 2 )---H I (T a ) 



and Hi (t) is given by 



Hj(r) : 



p+ 1 



1 




d~y(ty p I +1 {T,~y) 



Using the time ordered product, the operator U is given by the well known formula (see [12] p. 85): 



Hence we recover the Heisenberg formulation of interacting quantum field as a quantization of 
Butcher series. From another point of view we prove that formally the interacting quantum field 
given by the Heisenberg identity can be expressed as a Butcher series. 

1 Planar p— trees 

Let us introduce some definitions concerning planar trees. 

Definition 1.1 A planar tree is an oriented connected finite graph without loop together with an 
embedding into the plane; we suppose that the graph has a particular node that no edge points 
to; this node is called the root of the tree. The root is drawn at the bottom of the tree. 

Remark 1.1 

The set of planar trees differs from the set of rooted trees used by D. Krcimer and A. Connes in 
their work about renormalization [11], [7]. For instance the following planar trees are different 



although they represent the same rooted tree. All the results of this paper can be expressed using 
rooted trees instead of planar trees: we just have to adapt the definitions and add a symmetry 
factor in front of each term of the perturbative expansions. We prefer to work with planar trees 
because they permit to avoid these symmetry factors and to get simpler formulaes. 

Notation 1.1 

Let b be a planar tree 

1. The external vertices ofb are called leaves and the other vertices internal vertices, \b\ denotes 
the number of internal vertices of b. 

2. We denote by o the planar tree without internal vertex. 

3. A planar tree is a p-tree if each of its internal vertices has exactly p childrens. We denote 
by T(p) the set of p-trees. 

Example 1.1 

The planar trees of remark 1.1 are 2-trees and they satisfy \b\ = 2. 

Definition 1.2 Let (b\, . . . , b p ) <G T(p) p , then we denote by B + (b\, . . . , b p ) the p-tree obtained by 
connecting to a new root the roots of b\ and 62 and . . . and b p . 





Property 1.1 Let b e T be such that 6^0 then there exists a unique p-uplet (b\, . . . , b p ) G T(p) p 
such that b = B + (b\, . . . ,b p ). 



2 Perturbative classical field theory: Butcher series 

Let q be an integer such that q > d/2. Then it is well known (see e.g. [1]) that the Sobolev space 
H q (R d ) is an algebra, more precisely we have 



Proposition 2.1 Let f and g belong to H q (R d ) then the product fg still belongs to H q (R d ) and 
there exists a constant c q > which depends only on q and d such that 

ll/5ll_H"9(E d ) < c gll/lli/9(R ci )llsllff<!(R d ) 

From now on, for all q e N, H q denotes the Sobolev space H q ( 



(2.1) 



Let T > and (<p°, ip 1 ) € H q+1 x H q . Consider the following problem 

f <p<EC 1 ([0,T},H q )nC 2 ([0,T},H q - 1 ) 
(□ +m 2 )ip + Xip p = in H 9 ' 1 

<p(0,*) = <p° ; ^(O,.)^ 1 
We define recursively the family (</H^))beT(p) by setting 

0^)( t) t) := Sm ^V (fc) + cos(u k tyf(k ) (2.2) 

where for all ip e C°([0, T], R d ), ip(t, k ) denotes the spatial Fourier transform of ip. For all 
(61, . . . , b p ) e T(p)P wc set for all x G K" 

4>(B+(b 1 ,...,b p )){x):=- f G ret (x-y)cP(b 1 )(y)---<P(b p )(y) (2.3) 

where P + := {(t, a?) e R n |i > 0} and G ret denotes the retarded Green function of the Klein- 
Gordon operator i.e. 

G ret(z) := tA^CO / d d k S ' m{z0LJk) e^d n y (2.4) 

9 denotes the heavyside function (Q(t) = if t < and 1 otherwise). 
Remark 2.1 

We can easily check that <j>(o) is the solution of problem (2.1) with A = i.e. 

(□ + m 2 )cj>{o) = 

( o )( O,.) = ^;^(O,.) = ^ 

So we can consider <j>(p) as " the free field corresponding to the interacting field (f at time t = 0" . 
In other hand <j){B + (b\ . . . b p )) satisfies 

(□ + m 2 )0(B + (6 1 . . . bp)) = -</>(&!) • • • <t>{b p ) 

with zero Cauchy data on the hypersurface t = 0. 



Then we have the following result 

Theorem 2.1 For all T > and (/,/) € H q+1 x the family {4>{b)) bel{p) is well defined by 
(2.2) and (2.3) and Vfe e T(p), 0(6) beiongs to C 1 ([0, T], i? 9 ) n C 2 ([0, T],H q ^ 1 ). Moreover there 
exists a constant C > which depends only on m, d, p and q such that if 

C\X(1 + MT)\ + W^Wf- 1 <1 (2.5) 
(here M denotes the constant M := max(m, 1/m)) then the power series 

6GT(p) 

converges in the C 1 ([0, T},H q ) n C 2 ([0, T] 1 H q ^ 1 ) topology and the sum ip is a solution of problem 
(2.1). 

Remarks 2.1 

We have studied the equation (2.1) but our approach can be extended for analytic non linearity. 
In this case, all planar trees are involved and condition (2.5) must be adapted (see [9] for more 
details). 

Proof : (of theorem 2.1) 

The proof of theorem (2.1) is very simple. Suppose that the power series ip = J2beT(p) ^' & '<A(&) 
converges then a simple calculation show that we have ip(0, •) = </>(o)(0, •) + = / and §f (0, •) = 
^■{0, •) + = / and for all (h . . . b p ) e T(p) p 

(□ + m 2 )0(B+(6 1 , . . . , b p j) = -0(h) . . . 4>(b p ) 

But for all b G T(p), b ^ o there exists a unique p-uplet (b\...b p ) of p-tree such that b = 
B + (b\ . . .b p ) hence we get 

(□ + m 2 V=(n+m 2 )«Ko)+ ]T A^ 1 -- b ^(n + m 2 )<j)(B + (b 1 ,...,b p )) (2.6) 

(6i,...,6 p )eT(p)P 

=0-A \ lbll+ - +lbpl (f>(b 1 )---<j>(b p ) = -\ipP (2.7) 

(6i,...,6 p )€T(p)P 

Let us focus on the convergence of the power series. Let us show recursively that <p(b) belongs 
toC 1 ([0,T],ff'')nC 2 ([0,T], J ff«- 1 ) for all b e T(p) and 

11^(6)11 < ril + MT)) 1 ' 1 [M 2 (||^|| + ||^||)] |b|(p - 1)+1 (2.8) 

where M denotes the constant M := max(m, 1/m) > 1. 

Let us show (2.8) for b = o. The function 0(o) is given by (2.2). Then we can easily show that 
0(o) belong to C 1 ([0,T],i?")nC 2 ([0,T],i?9- 1 ) and since (to 2 + a 2 )/(l + a 2 ) < M 2 for all a e R 
we have 

||0(o)(t,»)lk« < M||/||h* + ||/||if* < M 2 (11/11 + 11/11) 
\\^-(t,»)\\m < y'Wm + M||/|| ff9+1 < M 2 (||/|| + ||/||) 

\\^l( tl .)\\ m -, < M||/|| H „ + M 2 ||/||^ +1 < M 2 (11/11 + 11/11) 
so ||0(o)j| :=max te[0 , T] (\\<j>(o)(t,.)\\ H «; ||^f%*)|k«; \\^(t,*)\\ Hq -i) satisfies (2.8). 



Suppose that (2.8) is satisfied for all b E T(p), |6| < N. Let b E T(p) such that |6| = N + 1 > 1 
then 3(bi...b p ) E T(p) p such that b = B + {b\ . . .b p ) and <j>(b) is defined by (2.3) hence, using 
proposition 2.1 we get 

IW&)(t,-)lk. < C^TH^oil • • • IWMII < (i + mt)^- 1 !!^)!! • • • 11^)11 



at 



-(t, < 05-^11^(61)11 • • • ||<KM|| < (i + MT)^u{b,)\\ ■ ■ ■ U{b p ) 

||^P(t, .)\\ Hq -r < <*~\\ + MT)\\4>{b,)\\ ■ ■ ■ H(b p )\\ 



hence we have ||<^(6)|| < (1 + MT)c^ _1 ||<^(6i)|| • • • U(b p )\\. Since (2.8) is satisfied for &i, ... , b p and 
using the fact that \B + (b\ . . . b p )\ = \b\\ H + \b p \ + 1 we finally get (2.8) for b. 

Finally it can be shown (see e.g. [16]) that the number of p-trees b such that \b\ = N is bounded 
by (joP/ip— l) p_1 ) Ar hence we see that if we have 

(1 + MT)|A| (11/11 + |^|)-< -JE^i£L (2.9) 

then the power series ^2 N ^1^11^(6)11 converges which completes the proof. ■ 
Remark 2.2 

Notice that the proof gives an explicit value for C. 

3 From Butcher series to perturbative quantum field theory 

We have seen that Butcher series provides an explicite and precise perturbative expansion of the 
solutions of interacting Klein-Gordon equation 

(□+m 2 V + A(^ = {E p ) 

So we get a perturbative classical field theory. 

On the other hand, physicists devollopped a perturbative quantum field theory (see e.g. [12], 
[10], [15]). The question we are interesting in is the following: is there a link between Butcher 
series and perturbative quantum field theory ? 

The answer is yes, we will show that the Heisenberg picture of interacting quantum field can be 
written as a quantized Butcher series. 

Let to E K. be a fixed time and x = (t, ~x) E M 1 ^. Then following M. Pcskin et D. Schroeder 
[12] p. 83, the free field on space-time is the self-adjoint operator 4>i(x) acting on Fock space (see 
[14], [13]) defined by 



4>i{x) = cj>i{t, ~x) := 



:=-!-/ -£L (3L k e- ik - x +ate ik - x ) 
(2nY J Rd ^T k \ k k ) 



x u =t-t 



Here k ■ x denotes the quantity k ■ x := LOkt — (k, ~~x) where (k, ~~x) denotes the usual scalar product 
on R d and ujk the quantity uju '■= (|fc| 2 + m 2 ) 1//2 . The operator and 3l\. denote the usual creation 
and annihilation operator (see e.g. [13]). 

We formally define the family (4>(b)(x)) be f^ of operators by the following 



Definition 3.1 <j>(o)(x) := cpi(x) and for all (h, ...,b p )e T{p) p 

0(B + (b u ...,b p ))(t,^):=- [ dy° f dl/G ret (x-y)$(b 1 )(y)--^(b p )(y) (3.1) 

Jto JM d 

Hence we can see the sum X)beT(p) A' fc '0(6)(:r) as the formal quantization of the Butcher series 
E 6eT ( P )A |61 ^). 

Remark 3.1 

• The definition 3.1 is formal because we do not care about the definition of the operator 
product which is well known to be ill defined. 

• We have the following commutation relation between the operators <j>i(x) and 4>i{y) 

A(x - y) := [M x), M y)] ~ (e"*^ - e*<-«<>) (3.2) 

(see [12] p. 28). Using the expression (2.4) of the retarded Green function G retl we get 

G ret (z) = -i9(T)A(z) (3.3) 
This simple remark leads to the following theorem: 
Theorem 3.1 for all t > t and ~~x <E R d we have 

(p(t,t) := J2 xlbl H b )( x ) = U\t) 4>i{t,lc) U(t) (3.4) 

beT(p) 



where U(t) denotes the operator 

U(t) := ]T (i\) a f dn f 1 dr 2 • • • f " 1 dTaH^H^Vi) ■ ■ ■ tf/(r Q ) 



(3.5) 



a >0 Jto Jto Jto 

and Hi(t) is given by 



y) 



Remarks 3.1 

• Using the time ordered product (see [12] p. 85 for a definition) the operator U is given by the 
well known formula 



U(t) =T jcxp ^A^" driJ/(r) 



• The right hand side of (3.4) is exactly the Heisenberg picture of interacting quantum field 
and this can be the beginning of the perturbative quantum field theory (see e.g. [12] p. 77-87). 

Proof : (of theorem 3.1) 

Let us introduce a new notation which will be very helpful for manipulations of iterated integrals 
such as (3.5). Let r £ N* and t < t- Then the symbol f> o dyi ■ ■ ■ dy r denotes the integrals 

/' d yi ---dy r := f dyl H dy° 2 ■ ■ ■ f'^ dy° p [ d d yl • • • f d d y* r 

Jto Jto Jto Jto JM d JR d 



Hence using this notation, expression (3.5) becomes 

u ® = E (^rr) Q ' ' • ■■4>T\y a ) (3.6) 

VP / J to 



So since 0/(x) is self-adjoint we get the following identity 

/ ' X \ " l 

[/t(^MMi)^ ^ E 

m>0 ^ ' (r,s)€N 
r-\-s=m 

(-If { dyi ■ ■ ■ dy r /' dzi ■ ■ ■ dz s ^ +1 (y r ) • • • ^ +1 (y 1 )0 / (^)^ +1 (^i) ■ ■ ■ <lF\z.) 

Jt Jt 

To complete the proof of theorem 3.1, it suffices to show that for all m € N 

nt ft 



P'/ ^ ,1cm "'to •'to 

# +1 (j/ r ) • • • ^ +1 ( 2/1 )0 / ( a; )^ +1 ^ 1 ) • • • 0? +1 (z s ) (3.7) 



beT(p) x ^ 7 (r,s) 

|6|=m r+s=m 



Let us prove (3.7) by induction on m e N. 

For m = identity (3.7) is obvious since the left hand side of (3.7) reduces to 4>i(x) and the only 
p-tree such that |&| = is given by b = o therefore the right hand side equals to 4>(o)(x) = <j>i(x). 
Let N e N and suppose that identity (3.7) is satisfied for all m < N. Then set 

ip N+ i(x) := Y &( b )( x ) 

beT(p) 
\b\=N+l 

Since N + 1 > 1, for all b € T(p) such that |6| = iV + 1 there is a unique (6i, . . . , b p ) G T(p) p such 
that & = B + (b\, . . . ,b p ). So the definition of 4>(b)(x) together with remark 3.1 lead to 

<Pn+i= Y 1 I d »° I d^A{x-y)^(b 1 ){y)--4(b p )(y) 
(6i,...,6 p )eT(p)f Jt ° jRd 
|6i|+-+|6 p |=JV 

which can be rewritten as 

vn + i= Y i f A y Q I ^A(x- y )( y kbi)(y))---( Y w P )(y))- 

qi-\ hq P =N |6i|=9i \b P \=q P 

Then using (3.7) for m < N, we see that <pn+i is given by the following expression 



N 

I 

I 



p+1 



E E E Mr 

( qi ,...,q p )&1 p ( ri ,si)eN 2 (r p ,s p )eN 2 
3iH hg P =A r ri+si=gi r p +s p =q p 



+-+r p 



f dy / dy^-y)/ dy« / dz« • • • / dy ^ / d^ p P((yf ),, fe , y) (3.8) 

Jto JM d Jto Jt Jt Jt 



Here we must explain the notations: 



• Let a, [3 be some integers such that a < (3, then we denote by y a ^ the (J3 — a + l)-uplet 

Va,P = (Va, ■ ■ , UP) € (R n ) /3_a!+1 . 

• IP((j/j' C ^)i,fe, y) denotes the product 

^ +1 (y^) ■ ■ ■ ^ + \y^)My) ■ ■ ■ ^SH^WSHyg) ■ ■ ■ <fi +1 (y[ 2) )] My) 

^ +1 (^ 2) ) • • • tf + \z?M + \^) ■ ■ ■ ^? +1 (4 3) )] My) ■ ■ ■ Mym +1 (z[ p) ) ■ ■ ■ 

Then we use the well known combinatorial lemma of quantum field theory (see e.g. [12]) which 
can be proved directly using iterated integrals 

Lemma 3.1 For all t > we have U{t)W{t) = Id i.e. for all m e N 
E (-If /' dy[ ...dy' r f d Vl ... dy s tf +1 ( yi ) ■ ■ ■ 

r-\-s=m 

where 5o, m = if m ^ and 1 otherwise. 
Using lemma 3.1 we see that expression (3.8) reduces to 



p I +1 (y'r)---r I +1 (y[) = s ,r, 



-i-rY E (-l) r /V/ d^A( 



x- y) 



(r.s)eN' 



i d Vl ---dyA &z 1 ---dz s ^ + \y r )---cj? I + \ Vl )My)€ + \zi)---€ +1 {zs) (3.9) 

J to ^io 

But since A(x — y) = [<fii(x),(f>i(y)] which commutes with <pi{z) (it is a c-number), we can replace 
A(x — y)4>i(y) p in (3.9) by its expression using 4>i{x) and <f>i(y) i.e. 

A(x - y)My) P = 0/(x)^ +1 (y) - Mv)M*)fl(v)- 
But we have <pi(x)<p p (y) = <p^(y)<pi(x) + pA(x — y)(f> I ^ 1 (y) therefore the last identity leads to 

A(x - y )Myf = M*Wi +1 {y) - 4> p+ \y)Mx) - pM* - y)Mv)- 

Inserting this last expression in (3.9) we finally get 

N+l 

E Mr, 

'to 



p+l 



t 

E (-l) r / dy° / dl?fdy 1 ---dy r fdz 1 ---dz s 



(r. 
r+s=N 



4> p+ \yr) ■ ■ ■ ^MMx^iyW^izi) ■ ■ • <tF\z 8 ) 

- ^ +1 (y I .)...^ +1 (y 1 )0? +1 (y)^(x)^ +1 (z 1 )...^ +1 (z s )] (3.10) 

Consider separately the terms of the right hand side of (3.10). The first term is given by (modulo 
a factor (i/(p + 1)) N+1 ) 



t 

E (-l) r / dy° / d^fdy 1 ---dy r fdz 1 ---dz 

s^ M 2 Jt JS. d Jto Jto 



(r,s)er 

r+s=N 



^ +1 {y r ) ■ ■ ■ 0? +1 ( 2/1 )0 / (^)^ +1 (y)^ +1 (^i) • • • 4> p+1 (z s ) (3.11) 



Notice that for all y € [0, t] x R d and a > 1 we have 

rt 

1 1 Jto Jy" ./!.-' Jt 



> dj/i • • • dy a =i dyi - • • dy a / dy° / dyl > dy 2 --- dy a 

Jto Jto Jy° JM. d Jto 



Hence using this last identity and performing the change of variable s <— s + 1, expression (3.11) 
leads to 

E h 1 ) 8 f d yi--- d Vrf dz 1 ---dz a Y(y r ...y 1 ,x,z 1 ...z s ) 

r>0;s>l 
r+s=7V+l 

^ (-If fdz* [ di? f dyl [ dyt 

n g>1 Jto JWi d Jz° JR d 



rvi r z i 

> dy 2 ---dy r -i dz 2 ■ ■ ■ dz a V(y r . . . yi, x, z\ . . . z s ) (3.12) 

Jto Jt„ 



!->l 

r+s=N+l 

,,o 



'to Jto 

where V(y r ■ ■ ■ yi, x, z\ . . . z s ) denotes the product 

V(y r . . .y u x, Zl . . . z s ) := 0? +1 (y r ) • • • ^ +1 ( yi )<M^ +1 (^i) ■ ■ ■ € + \zs) 

A similar computation (only replacing y by z and r by s) shows that the second term of the 
right hand side of (3.10) writes (modulo a factor —i I J ): 



r>l;s>0 
r+s=N+l 



r>l 
r+s=N+l 



E h 1 ) 3 f d yi--- d yrf dz 1 ---dz a V(y r ...y 1 ,x,z 1 ...z a ) 

E f d y" I d ^ f dz " I ^ 

>ls>1 Jto Jv, d Jy° Jw d 

> dy 2 ---dy r -± dz 2 ■ ■ ■ dz s Y{y r . . . yi, x, Zi . . . z s ) (3.13) 

Jto Jto 

Hence injecting (3.12) and (3.13) in (3.10), we finally find out that operator (^-j^-) N+1 fN+i is 
given by 

2 E ( _1 ) S f d yi--- d yrf dz 1 ---dz s V(y r ...y 1 ,x,z 1 ...z s ) 

r>l;s>l *° *° 

r+s=AT+l 

+ (-l)** 1 i dj/i • • • dyAT+xV^Ar+x . . . yi,a;)+ / d?i ■ ■ ■ dzjv+iV(ic, Zl . . . zjv+x) + A (3.14) 
where A denotes the operator 

- E (-W d ^ / arj [\z<{ C dy° 1+ C dz°) 

r>i~r>i ■ /Rd ■ /Rd V- 7 *" " / *o "M / 

> dy 2 ---dy r i> dz 2 ■ ■ ■ dZsV(y r . . . y 1 , x, z\ . . . z s ) (3.15) 

Jto Jto 



r+s=N+l 



But we have the following identity 

f&4 fdyU fdy* f*zl= fdz* fdy* 

Jt Jz\ Jto Jyl Jt a Jt 

So inserting this last identity in the expression (3.15) of A we get 

A=- ^2 " f cfyi • • • dy r > dzi---dz s V(y r ...y 1 ,x,z 1 ...z s ) 

r>l;s>l *° *° 

r+s=N + l 

hence we finally see that (3.14) leads to 

N + l rt 

l Pn+i= ( _1 ) S f dyi---dy r > dzi ■ ■ ■ dz s Y(y r ■ ■ ■ yi,x, z\ . . . z s ) 

r>0;s>0 ^ t( > 

r+s=N+l 

which is exactly (3.7) at order N + 1. ■ 
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